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Abstract 

In this work we compute, at the "one-loop-dressed" level, the nonperturbative contribution of the 
ghost loops to the self-energy of the gluon propagator, in the Landau gauge. This is accomplished 
within the PT-BFM formalism, which guarantees the gauge-invariance of the emerging answer. 
In particular, the contribution of the ghost-loops is automatically transverse, by virtue of the 
QED-like Ward identities satisfied in this framework. Using as nonperturbative input the available 
lattice data for the ghost dressing function, we show that the ghost contributions have a rather 
sizable effect on the overall shape of the gluon propagator, both for d = 3, 4. Then, by exploiting 
a recently introduced dynamical equation for the effective gluon mass, whose solutions depend 
crucially on the characteristics of the gluon propagator at intermediate energies, we show that 
if the ghost loops are removed from the gluon propagator then the gluon mass vanishes. These 
findings strongly suggest that, at least at the level of the Schwinger-Dyson equations, the effects of 
gluons and ghosts are inextricably connected, and must be combined suitably in order to reproduce 
the results obtained in the recent lattice simulations. 

PACS numbers: 12.38.Aw, 12.38.Lg, 14.70.Dj 
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I. INTRODUCTION 



Our understanding of the infrared (IR) properties of the fundamental Green's functions 
of Yang-Mills theories has improved considerably in the last few years, due to a va riety 
of parallel efforts in lattice simulations jl 12 1 , Schwinger-Dyson equations (SDEs) 13-20], 



functional methods 



2l|,|22), and algebraic techniques 23N26|. The majority of the aforemen- 



tioned studies have focused on the low-momentum behavior the gluon and ghost propagators, 
which can be directly or indirectly related to some of the most fundamental nonperturba- 
tive phenomena of QCD, such as quark confinement, dynamical mass generation, and chiral 
symmetry breaking. 

It is by now well-established that, in the Landau gauge, the lattice yields a gluon propaga- 
tor and a ghost dressing function that are finite in the IR (in d = 3, 4) 



14 



20j. Evidently, 



these lattice results furnish strong support to the idea of dynamical gluon mass generation 



through the well-known Schwinger mechanism 27H29I]. as proposed by Cornwall and oth- 



ers 



30 



31] . On the other hand, these important lattice findings have motivated the critical 



revision of the original Gribov-Zwanziger confinement scenario, leading to the formu 
its "refined" version 



. In addition, the "ghost-dominance" picture of QCD [32|, [33J , whose 



at ion of 



theoretical cornerstone has been the existence of a divergent ( "IR-enhanced" ) ghost dressing 
function, is at odds with the above lattice results, and, at least in this strict formulation, 
has been practice^ ral ed out f m the La lldaU gaug e, and ford = 3,4) QflQ. 

This last statement, however, does not necessarily mean that the ghost has been relegated 
to a marginal role in the QCD dynamics. In fact, compelling evidence to the contrary 
has emerged from detailed studies of the gap equation that controls the breaking of chiral 
symmetry and the dynamical generation of a constituent quark mass 34|, |35(. Specifically, 
the proper inclusion of the corresponding ghost sector (essentially the ghost dressing function 
and the quark-ghost kernel) is crucial for obtaining a realistic symmetry breaking pattern, 
with quark masses in the phenomenologically relevant range. The main lesson that can 
be drawn from the above studies is that even a finite (i.e., "non-enhanced") ghost sector 



may have a strong numerical impact, at least in the framework of the SD 



nontrivially the realization of various underlying dynamical mechanisms 35(. In fact, for 



Is, and affect 



the concrete case of the quark gap equation, the ghost contributions provide the necessary 
enhancement to the kernel of the gap equation precisely in the range of momenta around 



2 



1 GeV, which is the most relevant for obtaining the right type of quark mass solutions 36]. 

Given the importance of the ghost sector for the dynamical generation of a constituent 
quark mass, it is natural to ask whether a similar situation applies in the case of the dy- 
namical generation of an effective gluon mass. The main purpose of the present article is to 
address in detail this important question. 

This problem is technically rather subtle, and hinges on the ability to treat self- 
consistently various field theoretic ingredients. To that end, we will employ the general 
formalism based on the pinch technique (PT) 15|, |30|, |37H39| and the background field 



method (BFM) [40], which is particularly suited for dealing precisely with this type of prob 



lem. Specifically, the truncation scheme based on the PT-BFM formalism 41 



431 allows for 



subtraction of the ghost contributions to the gluon self-energy in a physically meaningful way 
(i.e., without introducing gauge artifacts). Indeed, in the conventional SDE formulation, 
any attempt to isolate the ghost contributions is bound to interfere with the transversality of 
the resulting gluon self-energy; this can be seen already at the one-loop level, where only the 
sum of the gluon and ghost diagrams (but not their individual contributions) is transverse. 
Instead, as was first pointed out in the classic paper by Abbott [40;], the calculation of the 
same diagrams using the BFM Feynman rules gives rise to two transverse contributions. 
This crucial property persists unaltered at the level of the SDE for the gluon self-energy: 
the SDE is composed by concrete subsets of "one-" and "two-loop dressed" diagrams, which 
are separately transverse, e.g., g^ITjf (g) = 0, where H% u (q) = (a 3 )^ + (a^)^ - see Fig. [TJ 
Therefore, one can study the individual contribution of the different blocks [in this case 
Il^(g)] to the full gluon self-energy, without compromising the transversality of the answer 
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41 



(these points have been addressed in great detail in 

Given that within the PT-BFM framework the ghost contributions to the gluon self- 
energy may be disentangled gauge invariantly, the next step will be to compute this particular 
contribution nonperturbatively, and then subtract it out from the full gluon propagator 
obtained from the lattice. The basic operating assumption underlying this analysis is that 
the gluon propagator found on the lattice coincides with that obtained from the solution of 
the full SDE series. Then, instead of solving the SDE series without the loops contained 
in U^ u (q) to determine the resulting gluon propagator (technically an impossible task at 
the moment), we compute nonperturbatively only the contribution of U^ u (q) and subtract 
it from the gluon propagator obtained from the lattice. The nonperturbative computation 
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FIG. 1: The SDE corresponding to the PT-BFM gluon self-energy H°^,(q). The graphs inside each 



box form a gauge invariant subgroup, furnishing an individually transverse contribution. White 
(black) circles denote full propagators (vertices). 

of the aforementioned ghost contribution H£ u (q) [graphs (03) and (04) in Fig. [Q proceeds 
through the following main steps. 

(i) We introduce a suitable Ansatz for the full (background) gluon-ghost vertex, T^, [the 
black circle in graph (03)] which satisfies automatically (i.e., by construction) the 
all-order Ward identity given in Eq. ( 12. 6p . This is an indispensable requirement for 
maintaining the gauge invariance (transversality) of the answer. The obtained from 
this procedure [given in Eq. (13.ip ] is expressed entirely in terms of the ghost propagator. 
As a result, the only quantity appearing finally inside the graphs (a 3 ) and (a 4 ) is the 
ghost propagator (or its dressing function). 



(ii) We invoke the so-called "seagull identity", given in Eq. (13 .4p . which enforces the can- 
cellation of all sorts of seagull-type contributions, leading to the absence of quadratic 
divergences |44j. Specifically, by means of this identity the purely seagull contribution 
of graph (a 4 ) cancels in its entirety against a term obtained from (03). 

(Hi) The remaining expression for H^ u (q) is renormalized subtractively, according to the 
rules of the momentum-subtraction (MOM) prescription. 
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(iv) The renormalized expression for IT^(g) is then computed numerically, by substituting 
for the (infrared finite) ghost dressing 
available lattice data for this quantity 



unction, appearing inside the integrals, the 



(v) The latter contribution is subtracted from the entire gluon propagator obtained from 
the lattice, according to the formula given in Eq. (I2.17p . 

The results turn out to be rather striking (see the right panels of Figs. El 151 and ITTj) : the 
gluon propagator without Il^ u (q) is significantly different from the full one. In addition, the 
results suggest a strong dependence on the space-time dimensionality: the effect of removing 
ghosts becomes considerably more enhanced as the space-time dimensionality is lowered. 

At this point one can turn to the main question of this work, and study what would 
happen to the gluon mass if the ghost contributions, computed in the previous steps, were 
to be removed from the full gluon propagator obtained from the lattice. This question 
can be addressed in quantitative detail by means of the integral equation, derived recently 



in 45], which describes the evolution (i.e., momentum-dependence) of the dynamical gluon 
mass, m 2 (q 2 ). This particular equation, given in Eq. (15. ip . contains as its main ingredient the 
full gluon propagator A, which practically determines the form of its kernel. The detailed 
analysis of an approximate version of Eq. ( 15. ip carried out in j^j] reveals that the existence 
of physically acceptable solutions hinges crucially on the shape of the gluon propagator in 
the entire range of physical (Euclidean) momenta, and in particular on the precise behavior 
that A displays in the region between (1-5) GeV 2 . Specifically, in order for the gluon mass 
to be positive definite, the first derivative of the quantity q 2 A(q 2 ) (the "gluon dressing 
function") must furnish a sufficiently negative contribution in the aforementioned range of 
momenta. Note that, as was shown in [45[, the full A obtained from the lattice has indeed 
this particular property, giving rise [when inserted into Eq. (15. ip ] to a dynamically generated 
gluon mass with the expected characteristics. Evidently, the main effect of removing the 
ghost contributions contained in II^(g) from A(g 2 ) is to restrict significantly the negative 
area displayed by the (ghostless) A(g 2 ) (see the right panels of Fig. [12]) . a fact which, in 
turn, leads to the vanishing of the gluon mass, i.e., the homogeneous Eq. (15. ip can only 
admit the trivial solution m 2 (q 2 ) = 0. 

Interestingly enough, this result appears to be completely analogous to what happens 
in the case of chiral symmetry breaking, where failure to include ghost contributions into 
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the gap equation [the quark analogue of Eq. (15.11) ] prevents the dynamical generation of 



a constituent quark mass 



351 ] . Thus, in the picture of QCD emerging from this analysis, 



ghosts and gluons must be in a state of harmonious synergy in order for a mass gap to be 
produced, regardless of the nature of the fundamental particle in question (gluon or quark). 

The article is organized as follows. In Section [III we study the general properties of the 
ghost sector in the Landau gauge and explain how it is possible within the PT-BFM frame- 
work to disentangle gauge- invariant ly the (one-loop dressed) ghost contributions, H^ u (q), 
from the full gluon propagator. In Section IIHI we derive the non-perturbative expression 
that determines I\.^ v {q) solely in terms of the ghost propagator and the coupling constant. 
In Section |IV] we evaluate numerically the expressions for IT^(g) derived in the previous sec- 
tion, using as input the ghost dressing function obtained in recent lattice simulations jl, 5, 9]. 
Next, we determine how the removal of H% u (q) affects the overall shape of the resulting gluon 
propagator, for three different cases: d = 4 and N = 2,3, as well as d = 3 and N = 2, where 
d is the dimensionality of space-time and N is the number of colors [corresponding to the 
gauge group SU(N)]. In Section [V] we turn to the main question of the present work, and 
study in detail how the kernel of the dynamical integral equation governing the gluon mass 
gets modified after removing the aforementioned ghost contributions. Finally, our conclu- 
sions are presented in Section [V]] 



II. GAUGE INVARIANT SUBTRACTION OF GHOST LOOPS 

In this section we first derive the formula that will determine the residual gluon propagator 
obtained from the full gluon propagator after removing from the latter the "one-loop dressed" 
ghost contributions, given by diagrams (a^) and (04) in Fig. [TJ Then, we work out the 
nonperturbative expression that determines the aforementioned ghost contribution in terms 
of integral involving the ghost dressing function. 

The first important fact to recognize is the transversality of the ghost contributions to 
be removed. Specifically, denoting their sum by 

nf (g) = (a„r + (a 4 r , (2.1) 
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we have that 

(a 3 U = -9 2 C A ff^{k,q,-k-q)D(k)D(k + q)V v (k + q,-q,-k), 
Jk 



(OA 



2g 2 C A g^ [ D(k). (2.2) 
Jk 



In the equations above, D ab (q 2 ) = 5 ab D(q 2 ) denotes the full ghost propagator, defined in 
terms of the ghost dressing function F as 



D{q 2 ) = ^P, (2-3) 



while T M represents the three-particle vertex describing the interaction of the background 
gluon with a ghost and an antighost, with (all momenta entering) 

tT cb ^ 5C (r, q ,p)=gr b f,(r, q , P y, f <?> (r, </, p) = (r-p) M . (2.4) 

Finally, Ca is the Casimir eigenvalue of the adjoint representation [Ca = N for SU(N)], 
and we have introduced the d- dimensional integral measure (in dimensional regularization) 
according to 

^ t d d A;, (2.5) 



k 



' (27T) 

with n the 't Hooft mass, and e = 4 — d. 

Then, by virtue of the PT-BFM Ward identity 



q »T fl (r } q,p) = D- 1 (r)-D- 1 (p), (2.6) 



it is immediate to establish the transversality of H% v (q), namely 41 j 



qJX? (q) = 0. (2.7) 

Let us now denote by H^ u {q) the sum of the remaining subsets of diagrams in Fig.HJ i-e., both 
the gluon one- and two-loop dressed diagrams, as well as two-loop dressed ghost diagrams, 

10 

nf(g) = J>*r- ( 2 - 8 ) 

i=l 

Again, due to the special Ward identities satisfied by the PT-BFM vertices, U% u (q) is 
also transverse, and, of course, so is the full self-energy given simply by 

nr(q) = n^( g ) + nf ( q ). (2.9) 




FIG. 2: Definitions and conventions of the auxiliary functions A and H. Gray blobs denote 1-PI 
kernels (with respect to vertical cuts). 

The SDE for the full gluon propagator in the Landau gauge of the PT-BFM scheme 
assumes then the form 



A- l (q 2 )P^(q) = 2 v ^ 1 r ^' ' c V ^ J , (2.10) 

where the gluon propagator A^(g) is defined as (we suppress color indices) 



A^q) = -zA(g 2 )P^(g); P, u (q) = g, v - (2.11) 

The function G appearing in (12.101) is the form factor associated with in the Lorentz 
decomposition of the auxiliary two-point function A, given by 

KM = ~W 2 C A I \l{k)D{q-k)H va {-q,q-k,k) 
Jk 

= 9»M<l 2 ) + qj ^Ltf). (2.12) 

This latter function, together with the auxiliary function H, are diagrammatically rep- 
resented in Fig. |2j also notice that H is related to the (conventional) gluon-ghost vertex by 
the identity 

p»H^{p,r,q) + r^r,q,p) = 0, (2.13) 
and that, in the (background) Landau gauge, the following all order relation holds [dfj 

F- l (q 2 ) = l + G(q 2 ) + L(q 2 ). (2.14) 

Now, let us return to Eq. (12.101) . and define in a completely analogous way the quantity 
A r (q 2 ), given by 
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Evidently, A r represents the propagator obtained by subtracting out (gauge invariantly) 
from the full propagator A the one-loop dressed ghost contributions. Then, taking the trace 
of both Eqs. (l2~T0j) and (j235j) . defining the trace of Ilf (?) as 

II c (g 2 ) ee n^(g), (2.16) 

and solving for A r , we arrive at 

I (d - 1) [1 + G( 9 ' 2 )] 2 / ' 

which represents our master formula. 

In order to obtain the behavior of the propagator A r (q 2 ) from Eq. ( 12.17ft we will (i) iden- 
tify the full gluon propagator A(g 2 ) with that obtained from the lattice, and (ii) determine 
nonperturbatively the quantity Il c from Eqs. ( 12. 2ft and (12.161) . and evaluate it numerically 
using as input the lattice results for the ghost dressing function F(q 2 ). These points will be 
the subject of the next two sections. 

III. THE NONPERTURBATIVE EXPRESSION FOR U c (q 2 ). 

To accomplish step (ii) above, we first need to introduce an Ansatz for the fully-dressed 
ghost vertex T^, appearing in graph a 3 of Eq. ( 12. 2p . which satisfies the crucial Ward identity 



of Eq. (12. 6p (this general procedure is known as the "gauge-technique 48]). The required 



Ansatz is easily constructed from that derived in 



49J for the case of scalar QED case, 



requiring the absence of kinematic or dynamical singularities. It reads 

r M (r,?,p) = [D^{ P 2 )-D-\r 2 )\ , (3.1) 

and evidently satisfies Eq. H 2 . 6 f) when contracted with q^. Obviously, the procedure of recon- 
structing the vertex by "solving" its Ward identity (known in general as "gauge technique" ) 



leaves the transverse (automatically conserved) part of the vertex undetermined [48l . 150 ]. 
In this case this term has the form A(r, q) [(r ■ q)p^ — (p ■ q)^^]. This particular term van- 
ishes as q — > 0, provided that the form factor A(r, q) does not diverge too strongly in that 
limit, which we will assume in what follows. Under this assumption, the transverse part 
of the vertex is subleading in the IR. On the other hand, its omission is known to affect 
the renormalization properties of the resulting SDE, a fact that forces one to renormalize 
subtractively instead of multiplicatively (see below). 



Substituting (13. ip in the first equation of (12. 2p and taking the trace, it is relatively 
straightforward to obtain the result 

n c (g 2 ) = g 2 C A [4T(g) - g 2 R(q)] , (3.2) 

where 

~D(k + q) -D(k) 



R(q) 



(k + q) 2 -k 2 



2 D(k + q)-D(k) d 



To further evaluate IT c (g 2 ), we must invoke the so-called "seagull-identity" [44J , 



valid in dimensional regularization, which enforces the cancellations of all seagull-type of 
divergences. This identity guarantees the (nonperturbative) masslessness of the photon in 
scalar QED [by setting f(k 2 ) -> V(k 2 ), where V(k 2 ) is the full massive scalar propagator], 
as well as the absence of quadratic divergences from the SDE determining the dynamical 



gluon mass [by equivalently setting f{k 2 ) — > A(k 2 )] [44j. 

For the case at hand, what we want to guarantee is that I1 C (0) = 0; this must be indeed so, 
because the ghost-loop giving rise to IT c (g 2 ) has no direct knowledge of the mass generating 
mechanism, namely the fact that A _1 (0) = m 2 (0). The easiest way to appreciate this is 
by recalling that the mechanism responsible for endowing the gluon with a dynamical mass 
relies on the presence of massless poles in the nonperturbative tree-gluon [the black circle 
in graph (ai) of Fig. dj, whereas the ghost vertex has the usual structure [note the absence 



of poles in the Ansatz of Eq. (13.11) ] [451 ] . 



Evidently, in the limit q — > 0, the term q 2 R(q) vanishes, and so does T(q), since 



= 0, (3.5) 

where in the last step we have employed Eq. (13 .4p . with f(k 2 ) — > D(k 2 ). 

In addition, note that the perturbative (one-loop) version of the terms R(q) and T(q), 
obtained from Eq. f]3. 31) by setting D{k 2 ) = 1/k 2 , is given by 

1 



RW(q) 

T«(g) = (d/2-1) I 4 (3.6) 



k 2 (k + q) 2 ' 
1 

Ik k 2 
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Evidently, due the dimensional regularization result J k k~ 2 = , we have that T^\q) = 0, 
and in the limit q — > 0, q 2 R^\q) vanishes (in d = 3,4). 

It is clear that, when d = 4, R(q) is ultraviolet divergent, and must be properly renormal- 
ized, by introducing in the original Lagrangian the appropriate counterterm or wave-function 
renormalization (the need to renormalize is seen explicitly already at the level of R^(q), 
which diverges logarithmically). The (nonperturbative) renormalization of U c (q 2 ) that we 
will employ proceeds as follows. First of all, as happens almost exclusively at the level 
of SDEs, the renormalization must be carried out subtractively instead of multiplicatively. 
The main reason for that is the mishandling of overlapping divergences due to the ambigu- 
ity inherent in the gauge-technique construction of the vertex, related with the unspecified 
transverse part |51| . 

The (subtractive) renormalization must be carried out at the level of (I2.10p . Specifically 
(setting directly d — 4), 

[1+G( 9 2 )] 

where the renormalization constant Za is fixed in the MOM scheme through the condition 
A _1 (/i 2 ) = fi 2 . This condition, when applied at the level of Eq. (13. 7ft . allows one to express 
Za as 

z A = [i + G(/i 2 )] 2 - ^ [rL» + n c ( M )] . (3.8) 



Now, as is well-known 



43, 



561 ]. the validity of the BRST-driven relation (I2.14p before and 
after renormalization prevents G(fi 2 ) from vanishing when, according to the MOM prescrip- 
tion, F(fi 2 ) = 1; instead, we must impose that G(fi 2 ) = — L(/i 2 ). However, given that L(x) is 
considerably smaller than G(x) in the entire range of momenta, we can use the approxima- 
tion 1 + G(fi 2 ) ~ F _1 (/i 2 ) = 1, without introducing an appreciable numerical error. Thus, 
we obtain the following approximate equation for Za 

Z A = l-^[U r (fi) + U c (fi)]. (3.9) 

Finally, substituting Eq. (I3.9P into Eq. (13. 7p . and defining (in a natural way) the renormalized 

A- X (g 2 ) as 

A r [q )- [1 + G(g 2 )] 2 ' [6A[)) 
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the renormalized version of the master formula (12. IT)) will read 

A-Ua 2 ) - A-Ua 2 ) - - [Hc(g) ~ (3 11) 

Evidently ( I3TTT]) is obtained from fl2TT7j) by replacing A^g 2 ) A^g 2 ) ("i?" for "renor- 
malized"), and II c (g) — > IT c ^(g), where 



U CiR (q) = U c (q) - (g 2 //i 2 )n c (^). 



(3.12) 



As an elementary check, note that the application of the last formula at one loop yields 

n^O?) = -g 2 C A q 2 [RV{q)-RV{n)] 

Vln(g 2 /^ 2 ), (3.13) 



ig 2 C A 2 . 



16tt 2 



which is the standard one- loop result of the PT-BFM 33, |40|, renormalized in the MOM 
scheme. 

For the ensuing numerical treatment of R(q) and T(q) carried out in the next section, it 
is advantageous to have the crucial property T(0) = a priori built in, in order to avoid 
possible deviations due to minor numerical instabilities. To that end, we introduce the 
quantity T 



T(q) = T(q) - T(0) 



k 



D(k + q) - D(k) _ dD(k) 

dk^~ 



(3.14) 



(k + q) 2 - k 2 

which has the property of ensuring (by construction) that T(0) = 0, while, at the same time, 
coinciding with the original T for all momenta q. 

In addition, it is convenient to re-express R(q) and T(q) in terms of the ghost dressing 
function. Using Eq. (I2.3p . after some elementary algebra, one obtains 



R(q) 

T(q) 



F(k) 



+ 



F(k + q)-F(k) 
k k 2 {k + q) 2 ' J k k 2 [{k + q) 2 -k 2 \ 
F(k + q)-F(k) dF(k)~ 



(3.15) 



(k + q) 2 - k 2 dk 2 

note that the angular integration of the first term in R can be carried out analytically for 
any value of the space-time dimension d. 

Finally, note that up until this point we have been working in Minkowski space. To 
make the transition to Euclidean space, we must employ the usual rules. Specifically, we set 
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1E-3 0.01 0.1 1 10 100 1E-3 0.01 0.1 1 10 100 



q 2 [GeV 2 ] q 2 [GeV 2 ] 

FIG. 3: Left panel: Lattice result for the SU(3) gluon propagator, A(q), in d = 4, renormalized 
at fi = 4.3 GeV. The continuous line represents the fit given by Eq. (|4.ip . Right panel: The 
SU(3) ghost dressing function, F(q 2 ), renormalized at the same point, fi = 4.3 GeV; the solid line 
corresponds to the fit given by Eq. ()4.3p . 

Ik = * Jfc E an d <?i = ~ < f i an d use that 

A E (<^) = -A(-gg); F E (ql) = F(-ql); G^ql) = G(-q 2 E ), (3.16) 
suppressing the subscript "E" in what follows. 

IV. NUMERICAL EVALUATION OF U c (q 2 ) AND A r (q 2 ). 

We will now proceed to perform the numerical analysis. Using the available lattice data 
on the ghost dressing function F, we evaluate the terms R and T given in Eq. (I3.15p . 
and combine them following the Eqs. (13. 2p and (13.121) to obtain the (renormalized) ghost 
contribution to the gluon self-energy Il c (of course, all relevant formulas must be properly 
"euclideanized" ) . Finally, we construct A r using (I2.17P and the lattice results available for 
the gluon propagator A. This exercise is carried out for three different cases: d = 4 and 
iV = 2, 3, as well as d = 3 and N = 2. 
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A. The case with d = 4, N = 3 



In Fig. [3] we show the lattice results for the four- dimensional SU(3) gluon propagator 
A(g 2 ) (left panel), and the corresponding ghost dressing function F(q 2 ) (right panel), ob- 
tained from [9], and renormalized at \x = 4.3 GeV. 



As has been discussed in detail in the literature 



35 



45 



52J, both sets of data can be 
accurately fitted in terms of IR-finite quantities. More specifically, for the case of A(g 2 ), we 



have proposed a fit of the form 



52] 



A-\ q 2 ) = M 2 (q 2 ) + q 2 



\ , ljjgAg? . fq 2 + PiM 2 (q 2 ' 
1 H zr—^ In 



96tt 2 



where 



M 2 {q 2 



m 



q 2 + p 2 m\ 



(4.1) 



(4.2) 



Notice that in the above expression, the finiteness of A _1 (g 2 ) is assured by the presence 
of the function M 2 (q 2 ), which forces the value of A _1 (0) = M 2 (0) = m\j ' p 2 . The continuous 
line on the left panel of Fig. [3] corresponds our best fit, which can be reproduced setting 
m = 520 MeV, g\ = 5.68, p 1 = 8.55 and p 2 = 1.91. 

The SU(3) lattice data for F(q 2 ), shown in the right panel of Fig. El will be fitted by the 
following expression 



F-\q 2 



1 | 9C ^j ln fQ 2 + PsM 2 (q 2 ' 



4 48tt 2 



M 2 {q 2 



ml 



q A + p 2 m^ 



2 - 



(4.3) 



with the parameters given by mo = 520 MeV, g\ = 8.65, p 2 = 0.68 and p^ = 0.25 . Notice 
that the M(q 2 ) has the same power-law running as the one reported in Eq. (14.21) 53M55| . 

It is interesting to notice that the aforementioned fits share the following important prop- 
erties: (i) they connect smoothly the IR and UV regions by means of a unique expression; 
(ii) their finiteness is associated with the presence of the parameter M in the argument 
of the perturbative (renormalization group) logarithm, which it is responsible for taming 
the Landau pole and for doing the logarithm saturates at a finite value 52j; and (in) for 
large values of q 2 , Eqs. (14. ip and (14. 3p reproduce their respective one-loop expressions in the 
Landau gauge. 

The only missing ingredient for the actual nonperturbative determination of Il c , and 
therefore A r , is the value of a s = g 2 /47r. Instead of choosing a single value for a s , we will 
establish a certain physically motivated range of values, which will furnish a more repre- 
sentative picture of the numerical impact of the ghost corrections on the gluon propagator. 
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Ghost dressing function 

EB V=64* and p=5.7 
a V=80 , and p=5.7 
SDE solution with a =0.29 




FIG. 4: The solution of the SDE (|4.4jl that best matches the ghost dressing function data is 
obtained for a s = 0.29. 

The lower value for a s will be fixed simply by resorting to the the 4-loop (perturbative) 
calculations in the MOM scheme [l. 9j , and extracting the value of a s that corresponds to 
the subtraction point of \x = 4.3 GeV, used to renormalize the lattice data. The value so 
obtained is a s = 0.2. 

In order to establish a reasonable upper bound, in a consistent way, we resort to the 



methodology employed in [561], which makes use of the standard SDE for the ghost dressing 



function, given by (Euclidean space) 
F-\q 2 ) = l + g 2 C A 







f k (k + q)* 


k 2 q 2 



A(k)F(k + q), 



(4.4) 



derived in the Landau gauge, and under the assumption that the full ghost-gluon vertex is 
approximated by its tree-level value jssl, I57]]. In this integral equation one substitutes for 
A(k) the fit given in Eq. ( 14. 1ft . and solves it numerically for the unknown function F(q 2 ); 
evidently, for each value of a s we obtain a different solution for F(q 2 ). The correct value of 
a s is then determined as the one for which the corresponding (renormalized) solution best 
matches the lattice results (see Fig. 0|); for \i = 4.3 GeV we obtain a s = 0.29, showing that 
the perturbative MOM value (a s = 0.2) is 30% lower. 

The results obtained for the renormalized R and T, after substituting into the corre- 
sponding formulas our best fit for F, given by Eq. ( 14. 3ft . are shown on the left panel of 
Fig. [51 together with the combination q 2 R — AT, which appears on the rhs of Eq. (13.21) . 
It is clear that the contribution of the term AT is rather negligible; in a way this is to be 
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q 2 [GeV z ] q 2 [GeV 2 ] 



FIG. 5: Left panel: Numerical evaluation of the ghost contribution H c (q) to the gluon propagator 
using as input our best fit for the d = 4, TV = 3 ghost dressing lattice data. Right panel: The 
removal of the one-loop dressed ghost contribution from the (lattice) gluon propagator results in a 
diminished "swelling" in the momentum region below 1 GeV 2 . 

expected, given that this term vanishes identically in perturbation theory (for all values of 
q), and vanishes nonperturbatively at the origin [viz. Eqs. (13.61) and (13. 5ft . respectively]. 

Next, we use these results to construct Il c , given in Eq. (I3.12p . and finally A r , expressed 
by Eq. (13. lip (Fig. [5] right panel), using both values of a s , namely a s = 0.29 (SDE, red 
dotted line) and ot s = 0.20 (4-loop MOM, blue dashed-dotted line). 

We then see that the net effect of removing the ghost contribution is to suppress signifi- 
cantly the support of the gluon propagator in the region below 1 GeV 2 . Higher values of a s 
increase the impact of the ghost contributions, but only slightly, as can be seen on the right 
panel of Fig. |5j As we will see in the next section, this "deflating" of the gluon propagator 
in the intermediate region of momenta, produced by the removal of the ghost contributions, 
has far-reaching consequences on the generation of a dynamical gluon mass. 

B. The case with d = 4, N = 2 

It turns out that, changing the gauge group to SU(2) does not significantly alter the 
characteristic qualitative behavior found in the SU (3) case. Specifically, in Fig. [6] we show 
the gluon propagator (left panel), and the ghost dressing function (right panel), obtained 
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q 2 [GeV 2 ] q 2 [GeV 2 ] 



FIG. 6: Left panel: Lattice results for the SU(2) gluon propagator in d = 4, renormalized at fj, = 2.2 
GeV. The continuous line represents our best fit to the data obtained from Eq. (|4.ip . Right panel: 
The SU{2) ghost dressing function F(q 2 ), renormalized at the same point, [i = 2.2 GeV; the solid 
line corresponds to the best fit given by Eq. (14.31) . 

from [7f and renormalized at /i = 2.2 GeV. 

As in the SU(3) case, the gluon and ghost data can be accurately fitted by the expres- 
sions (14. ip and (14. 3p . where now Ca = 2 and the fitting parameters are m = 865 MeV, 
g\ = 10.80, pi = 1.96 and, p 2 = 2.68 (gluon) and g\ = 15.03, m = 523 MeV p 3 = 0.215 and 
p 4 = 0.781 (ghost). 

The coupling a s can be also fixed using the same procedure described in the previous 
subsection (Fig. [7]); the value obtained from the SD solution that best matches the lattice 
data is in this case is a s = 0.99. 

On the left panel of Fig. [HJ we show the resulting curves for R and T obtained through 
our best fit for F given by Eq. (14. 3P . Then, using Eqs. (13 .21) and (13. IIP we combine the 
previous results to get the SU (2) ghost self-energy Il c and, finally, A r (right panel of the 
same figure). We use again two values for a s namely the one obtained through the solution 
of the ghost SDE (a s = 0.99) and a 30% lower one (a s = 0.70). Evidently, the SU(2) results 
do not differ qualitatively from those of the SU(3) case: a lower value for a s suppresses 
the ghost contribution to the gluon propagator, and the removal of the ghost gives rise to a 
lower curve in the region below 1 GeV 2 . 
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FIG. 7: The solution of the SDE (|4.4jl that best matches the ghost dressing function data is 
obtained for a s = 0.99. 




> 

CD 

CD 



V=128* and [1=2.2 
A(q) 

A r with ^=0.99 
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FIG. 8: Left panel: Numerical evaluation of the ghost contribution II C to the gluon propagator 
using as input our best fit (14. 3p for the d = 4, N = 2 ghost dressing lattice data. Right panel: The 
removal of the one-loop dressed ghost contribution from the (lattice) gluon propagator causes, as 
in the SU(S) considerable suppression in the momentum region below 1 GeV . 



C. The case with d = 3, JV = 2 

Let us start, as in the previous cases, by showing in Fig. 0the lattice results for 
the three-dimensional gluon propagator A(q) (left panel) and the ghost dressing function 
F(q) (right panel). Notice that, in Fig. [91 the lattice data for A(q) presented in Ref. jf], 
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FIG. 9: Left panel: Lattice results for the SU(2) gluon propagator in d = 3. The continuous line 
represents our best fit to the data obtained from Eq. (|4,5p . Right panel: Lattice data for the SU(2) 
ghost dressing function F(q) in 3 dimensions; the solid line corresponds to the best fit given by 
Eq. flM]). 

were appropriately rescaled, following the procedure explained in detail in [yj], to match 
correctly the perturbative tail. Both A(q) and F(q) saturate in the deep IR region, and can 
therefore be fitted by means of IR finite expressions. 

In the case of the gluon propagator, an accurate fit is giving by 

1 



A(q) = Aexp [-(<?- g ) 2 H + 



(4.5) 



a + bq + cq 2 ' 

where the fitting parameters are A = 0.49, q = 0.11, w = 0.37, a = 0.43, b = —0.85, and 
c = 1.143. For the ghost dressing function, we use the following piecewise interpolator 

1 



F(q) 



for q < 3 



a + bq + cq 2 ' 
1 + ~~ — t; ) for q 2 > 3 



(4.6) 



eq + q 2 

with fitting coefficients a = 0.19, b = 0.61, c = —0.14, d = 0.63 and e = 0.26 obtained by 
requiring the function to be continuous at q 2 = 3. 

The contribution of R and T of Eq. (13. 15ft can be then evaluated using the above fit, and 
the results of this calculation are shown in the left panel of Fig. [TTJ Since d = 3 Yang-Mills 
is a super-renormalizable theory, all aforementioned quantities are directly UV finite, and 
do not need to undergo renormalization. 
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FIG. 10: The solution of the SDE (14, 4p that best matches the ghost dressing function data in d = 3 
is obtained for g = 1.208. 

The next step is to determine the value of the coupling constant g (which, in d = 3, 
has dimensions of m 1//2 ) entering in the formulas for II C and A r , given by Eqs. (13. 2p and 
(I2.17p . respectively. The procedure followed is the same as before, i.e. we will employ 
the three-dimensional ghost SDE, solve it for various values of g, and choose the one that 
best reproduces the lattice data for F. The most favorable case is shown in Fig. [TOj where 
the solution for F(q) obtained from the SDE with g = 1.208 [in the same arbitrary mass 
units used in the plots of Fig. [9] (red line) is compared with the lattice results for the same 
quantity. 

Next, substituting the results presented on the left panel of Fig. [11] into Eqs. (I3.2p and 
(I2.17p . and using g = 1.208, we compute II C and A r . On the left panel of Fig.[TTl we compare 
the residual propagator A r (blue dashed line) with the full propagator A(g). Clearly, the 
effect in the tridimensional case is even more pronounced: the ghost contribution completely 
dominates over the rest, determining to a large extent the overall shape and structure of the 
propagator. 

V. NO GLUON MASS WITHOUT GHOST LOOPS 

In the previous section we have studied how the subtraction of the ghost contributions 
affects the profile of the gluon propagator. However, as we will now show, the effects goes 
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FIG. 11: Left panel: Numerical evaluation of the ghost contribution II C to the gluon propagator 
using as input our best fit for the d = 3, N = 2 ghost dressing lattice data. Right panel: The 
result of removing the one- loop dressed ghost contribution from the gluon propagator in d = 3. 
The effect is much more dramatic than in the d = 4 case, since all the structure is determined by 
the ghost contribution, while A r has the sole (but crucial!) role of rendering the propagator finite 
at q = 0. 

way beyond a simple change in the overall propagator shape, modifying its salient qualitative 
characteristics, and in particular the generation of a dynamical gluon mass. 

To establish this, we start from the dynamical equation describing the effective gluon 



mass, recently derived in 

2g 2 C A 



m 2 (q 2 ^ 



451 ]; it reads (Euclidean space) 

(k ■ q) 2 ] m 2 (k + q) — m 2 {k) 



k 



Q 



l + G(q 2 )Jk 
Taking the q — > limit, one then gets 

2g 2 C A d-l 



(k + q) 2 - k 2 



A(k)A{k + q). (5.1) 



m 2 (0) 



1 + G(0) d Jk 
4g 2 C A d-l 



k 2 [m 2 (k)}'A 2 (k) 
m 2 (k)A(k) [k 2 A(k)]' . 



(5.2) 



1 + 67(0) d Jk 

where in the last step we have used integration by parts. Introducing spherical coordinates 
(setting y = k 2 ) and the <i-dimensional integral measure [notice that in ( 15. ip there is no 
dependence on the d — 2 polar angles (fi] 

/ 1 71 2 

21 



I d9sm d - 2 e / dyy^ 
o Jo 



(5.3) 



Eq. (15. 2p finally becomes 



™ 2 (o) = - -Aj^tts r<W(v)£**(i/), (5.4) 



d(47r)fr(01 + G(O) J0 
with the kernel JC^n given by 

/C, ;J v(i/) = ^- 1 A(y)[yA(y)]'. (5.5) 

The dependence of JCd-N on iV (the number of colors) is implicit in the form of A(y) that 
must be employed in each case, i.e., A(y) = A N (y). The same is true for G(0) in (15 .4p . and, 
of course, Ca = N. 

Since the constant multiplying the integral is positive, the negative sign in front of 
Eq. (15 Ah tells us that the required physical constraint m 2 (0) > can be fulfilled if and 
only if the integral kernel Kd-,N (constructed solely out of the gluon propagator) displays a 
sufficiently deep and extended negative region at intermediate momenta 45| . 

In the left panels of Fig. [12] we plot the kernels ICd-N obtained from the lattice data for 
the cases d = 4, N = 3 (top row) and N = 2 (middle row), as well as d = 3 N = 2 (bottom 
row), considered in the previous section; they all posses the characteristic negative region 
that allows, at least in principle, the existence of solutions of Eq. (15.11) . furnishing a positive 
value for the condition (15. 2L We emphasize that, for the d = 4 and N = 3,2 cases such 



a solution has been explicitly found and studied in [45J. Notice the striking resemblance 
between the kernels obtained for the different cases. 

On the other hand, the situation changes substantially once the ghost loop is removed, in 
which case the kernels ICd-N must be constructed from A r (right panels of the same figure). 
For d = 4 one observes a shift towards higher gs of the zero crossing, and a correspondingly 
suppressed negative region; even though this is not sufficient to exclude per se the existence 
of a physical solution to the mass equation (15.11) . a thorough study of the approximate 
equation derived in |45| reveals that no physical solution may be found. The d = 3 situation 
is even more obvious: the highly suppressed negative region present in this case cannot 
support solutions of (15. 2\\ with m 2 (0) > 0, thus leaving as the only possibility the trivial 
m 2 = solution. 

The main conclusion one can draw, therefore, is that the ghosts play a fundamental role 
in the mechanism of dynamical gluon mass generation, since the failure to properly include 
them results in the inability of the theory to generate dynamically a mass for the gluon. 
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FIG. 12: The kernel K.^n of Eq. (|5.5p constructed out of the lattice propagator A (left panels) and 
the ghost-less propagator A r (right panels) for the d = 4 N = 3 (top row), N = 2 (middle row) 
and d = 3 N = 2 (bottom row) cases. The insets show in each case the shape of the propagator 
used to evaluate the kernels. 

This, in turn, implies that what is displayed in the right panels of Figs. El El EH and [12] are 
not the gluon propagators one would actually obtain, assuming that one were actually able 
to perform this "experiment", e.g., remove the ghosts on the lattice. Indeed, according to 
our results, if ghosts were not included, m 2 = 0, and thus A would not saturate in the IR 
at all! 
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To understand what happens, let us concentrate on the d = 4, N = 3 case and imagine 
a simplified setting, where one can switch off adiabatically the ghosts, neglecting all other 
effects this operation would entail (we will come back to this point at the end of this section). 
This could be achieved by multiplying the self-energy Il c , appearing in Eq. (|2.10|) . by a 
parameter 7 £ [0,1], such that when 7 = 1 the full A of the right panel of Fig. is 
reproduced. Now, by slowly decreasing 7 (for fixed a) one would give rise to a set of 
intermediate A rj7 profiles, showing progressively less "swelling" in the q 2 < 1 [GeV 2 ] region, 
and ideally one would get A ri0 = A r . However before that will happen, there will exist 
a critical value 7 C for which the kernel K,^ constructed from A r 7c will fail to provide the 
required negative region that would ensure the positivity of m 2 (0), as calculated from the 
condition (15 .4p . At that point the theory will undergo a drastic change, showing a gluon 
propagator that does not saturate in the IR. Even though we cannot actually predict what 
such a propagator might behave like in the deep IR, it is likely that the typical singularity 
associated with the (perturbative) Landau pole (tamed by the presence of the mass) may 
reappear. 

Obviously in this analysis we are neglecting any type of back-reaction due to the changes 
in the gluon propagator: to be sure, any modification to the latter quantity would affect not 
only the ghost - since the gluon propagator appears in fact in the ghost SDE, see Eq. (I4.4p 
- but also the gluon mass, and therefore the IR saturation value - through Eqs. (15. ip and 
(15. 2p . While such effects might be numerically appreciable (changing, e.g., the critical value 
7 C ), we expect the qualitative description given above to persist. 

VI. CONCLUSIONS 

In this article we have presented a detailed study of the impact of the ghost sector on the 
overall form of the gluon propagator in a pure Yang-Mills theory, for different space-time 
dimensions (d = 3, 4) and SU(N) gauge groups (N = 2, 3). 

The key ingredients for performing this analysis have been basically two. To begin with, 
the PT-BFM framework allowed us to subtract out gauge-invariantly the "one-loop dressed" 
ghost diagrams from the SDE describing the full gluon propagator. Second, we have been 
able to express these ghost contributions as a simple integral involving the ghost dressing 
function only. This was achieved by employing a judicious Ansatz for the ghost-gluon vertex, 
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obtained by solving the corresponding Ward identity, and by resorting to the "seagull- 
identity" , in order to enforce certain crucial properties. The nonperturbative evaluation of 
the resulting expressions have been carried out numerically, using available lattice data as 
input for the ghost dressing function. Our results reveal that the ( "one-loop dressed" ) ghost 
diagrams furnish a sizable contribution to the gluon propagator in d = 4, and the dominant 
one in d = 3. 

The suppression of the gluon propagator induced by the removal of the ghost-loops has 
far-reaching consequences on the mechanism that endows gluons with a dynamical mass, 
associated with the observed IR-finiteness of the gluon propagator and the ghost-dressing 
function. Specifically, using a recently derived integral equation controlling the dynamics 
of the (momentum-dependent) gluon mass, we have demonstrated that when the reduced 
gluon propagators are used as inputs, the corresponding kernels are modified in such a way 
that no physical solutions may be found, thus failing to generate a mass gap for the pure 
Yang-Mills theory. Instead, as has been shown in [45J], the use of the full gluon propagator 
in the same equation generates a physically acceptable gluon mass. 

Once the results of the present work are combined with those of [35J for the chiral sym- 
metry breaking, a compelling picture of QCD emerges, where the generation of a dynamical 
mass for quarks and gluons requires the synergistic participation of all fields (physical and 
unphysical) of the theory. 
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